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ABSTRACT 

Using very high resolution (A/AA ~ 95 000) 4.7 /xm fundamental and 2.3 overtone ro- vibrational 
CO absorption spectra obtained with the CRIRES infrared spectrometer on the Very Large Telescope 
(VLT), we report detections of four CO isotopologues - C^^O, ^^CO, C^^O and the rare species, C^^O 
- in the circumstellar environment of two young protostars, VV CrA and Reipurth 50. We argue that 
the observed CO absorption lines probe a protoplanetary disk in VV CrA, and a protostellar envelope 
in Reipurth 50. All CO line profiles are spectrally resolved, permitting direct calculation of CO oxygen 
isotopologue ratios with 5-10% accuracy. The ro-vibrational level populations for all species can be 
reproduced by assuming that CO absorption arises in two temperature regimes. For both objects, 
12(-i^13q g^j.g ^Yi^ order of 100, nearly twice the expected interstellar medium (ISM) ratio. The 
derived oxygen abundance ratios for the VV CrA disk show a significant mass-independent deficit of 
C'^^O and C^^O relative to C^^O compared to ISM baseline abundances. The Reipurth 50 envelope 
shows no clear differences in oxygen CO isotopologue ratios compared with the local ISM. A mass- 
independent fractionation can be interpreted as being due to selective photodissociation of CO in the 
disk surface due to self-shielding. The deficits in C^^O and C^^O in the VV CrA protoplanetary disk 
are consistent with an analogous origin of the ^^O variability in the solar system by isotope selective 
photodissociation, confirmation of which may be obtained via study of additional sources. The higher 
fractionation observed for the VV CrA disk compared with the Reipurth 50 envelope is likely due to 
a combination of disk geometry, grain growth, and vertical mixing processes. [Abstract abridged] 

Subject headings: astrochemistry — circumstellar matter — ISM: abundances — planetary systems: 
protoplanetary disks — stars: individual (VV CrA, Reipurth 50) 
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1. INTRODUCTION 

The capability to compare isotope ratios of meteorites 
with those in protoplanetary systems presents exciting 
new possibilities for constraining early solar system his- 
tory. Oxygen is of particular value in this regard because 
it occurs in great abundance, comprising ~ 50% by mass 
of most rocky materials and is, with the [CO]/[H] ratio 
roughly lO"'* in the gas phase, a significant portion of 
the gas as well. With the advent of improved instru- 
mentation, astronomical observations of stable isotope 
ratios are becoming a potentially powerful new tool for 
probing the chemistry of young stellar objects (YSOs), 
making comparisons with the solar system possible. 

Mass-dependent isotope fractionation describes the re- 
lationship between the changes in the ratios of three or 
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more isotopes. In the case of oxygen, one considers the 
relationship between the fractionation factors a^^/^^ and 
^18/16^ where a^^/ie ^ {[^^0]I^^O])l{[^^0]/[^^0]). and 
^18/16 ^ ([i80]/i60])/(rO]/[i60])o. The subscript o 
refers to some initial condition. Because mass-dependent 
fractionation depends on the ratio of partition func- 
tions, it is well known that the relationship is a^^/^^ — 
(q,i8/i6)/3^ where (3 = (l/mie - I / mn) / {I / miQ - l/mis) 
at equihbrium, or ln(Mi6/Mi7)/ln(Mi6/Mi8) for kinetic 
processes, and where rrii is the atomic mass for atomic 
species i, and Mi can be an atomic, molecular, or reduced 
mass for isotopologue i depending on the process (Young 
et al. 2002). In practice, /3 ~ 0.51 (kinetic) to 0.53 
(equilibrium) for oxygen. Isotope variations that obey 
this relationship are deemed mass-dependent. Mass- 
independent variations are those that do not follow this 
relationship. 

Oxygen isotopes in rocks of the solar system exhibit an 
anomalous mass-independent distribution that has defied 
conclusive explanation since its discovery (Clayton et al. 
1973). This anomaly, in the form of a mass- independent 
correlation between [i'50]/[^*0] and [^^0]/[^'^0] among 
rocky bodies, is one of the most pronounced chemical fea- 
tures of the solar system. Here, mass independence refers 
to relative differences in [160]/[180] and [i60]/[i^0] 
among primitive rocky objects that are nearly identical, 
suggesting changes in [^^O] relative to both [^''O] and 
[^*0], rather than the expected mass-dependent trend in 
which relative changes in [-'^^0]/[^^0] are about half those 
in [160]/[180] on a log scale. 
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Explanations for the mass-independent distribution 
of oxygen isotope ratios in the solar system include 
inheritance of a nucleosynthesis signal from the in- 
terstellar medium (e.g., Galactic chemical evolution), 
symmetry-based mass-independent isotope fractionation 
effects such as those involving ozone in Earth's strato- 
sphere (Thiemens & Heidenreich 1983), and photochem- 
istry involving isotope-specific photodissociation of CO 
(Kitamura & Shimizu 1983; Thiemens & Heidenreich 
1983; Clayton 2002; Yurimoto & Kuramoto 2004; Lyons 
& Young 2005). Recent results from the Genesis mission 
establish that rocks in the solar system are depleted in 
^•^0 relative to the Sim (McKcegan et al. 2008). The 
depiction is extreme in aqucously altered primitive me- 
teoritical materials (Sakamoto et al., 2007), suggesting 
that [Hi80]/[Hi60] and [Hi^Ol/pi^O] were extremely 
high in the planet-forming region of the early solar sys- 
tem. Both of these observations lend support to the sug- 
gestion that self-shielding by CO oxygen isotopic species 
during photodissociation of CO is a likely explanation 
for mass-independent fractionation of oxygen isotopes in 
the solar system (Young ct al. 2008). 

Self-shielding by CO refers to the variable shielding of 
CO oxygen isotopologues from photodissociation by far 
ultraviolet (FUV) radiation in proportion to their abun- 
dances (column densities). Thus, C^''0 and C^^O will 
be more rapidly destroyed than the much more abim- 
dant C^®0. The oxygen liberated during this process will 
eventually end up in H2O, providing an explanation for 
the overabundance of Hj^O and Hj^O relative to Hj^O 
in the early solar system. Exchange of oxygen isotopes 
between newly formed rock and H2O gas provided a path 
for rocks to become enriched in ^*0 and ^^O relative to 
the original solar oxygen isotopic composition (Yurimoto 
k Kuramoto 2004; Lyons & Young 2005; Young 2007). 

The effect on the survival of oxygen isotopologues is 
greatest for column densities of CO in the range 10"'^^ 
to 10^* cm^^. Isotope selective photodissociation and 
self-shielding occur in molecular clouds (Bally & Langcr 
1982; van Dishoeck & Black 1988; Sheffer et al. 2002) 
but their occurrence in circumstellar disks has yet to be 
demonstrated conclusively. Models of CO photodissocia- 
tion in a disk suggest that isotope selectivity will be most 
efficient at the surfaces of the disk regardless of whether 
the FUV source is the central star itself or proximal O or 
B stars (Lyons & Young 2005; Young 2007). These mod- 
els predict that the outer regions of circumstellar disks 
should exhibit C^^O and C^*^0 deficits relative to C^^^O 
as a consequence of CO self-shielding, providing a test of 
the models. The precise degree of isotope selection de- 
pends on the amount of radial and vertical mixing in the 
disk, as well as on time. Models suggest C-'^'^O and C-'^^O 
deficits relative to C^^'O of tens of percent on time scales 
of 10^ to 10^ years (Young 2007). This is more than suf- 
ficient to explain the anomalous distribution of oxygen 
isotopes in the solar system. An alternative model (Yu- 
rimoto & Kuramoto 2002; Lee et al. 2008) suggests that 
signatures of CO self- shielding in the disk may have been 
inherited from the surrounding cloud material, where the 
excess ^^O and ^^0 has subsequently been incorporated 
into water ice and transported to disks; in this case, one 
expects anomalous deficits in C^''0 and C^^O in the disk 
and in the envelope surrounding the disk, regardless of 
stage of evolution. 



In this work we present column densities for the CO 
isotopologues 12C160 (Ci^O), ^^C^^O (C^^O), ^^C^^O 
(C^^O), and "Ci^O (^^CO) in two YSOs that comprise 
a first test of the importance of CO self-shielding in the 
early solar system. We emphasize that quantifying both 
[i^0]/[iS0] and [^'^0]/["0] is particularly important in 
this study because only the combination of the two ratios 
can distinguish conclusively between photochemical self- 
shielding and mass-dependent fractionation. 

Infrared absorption spectroscopy has been used pre- 
viously to measure the ratios of the ^^CO, ^^CO and 
C^^O isotopologues in YSOs. Infrared observations 
of one circumstellar gaseous disk provide a hint of 
C^^O overabundance relative to C^^O: Brittain et al. 
(2005) presented high-resolution infrared spectra of the 
embedded, low-mass pre-main sequence star HL Tau, 
showing a ratio of column densities for C-'^O and 
Ci«0, iV(Ci60)/iV(Ci«0), of 800±200 {2a). Typical 
ji6o]/[i8o] in the local interstellar medium (ISM) is 
557 ±30, (Wilson 1999), similar to the solar system value 
of 499. The apparent overabundance of C^^O in the HL 
Tau disk could be the result of isotope-selective photodis- 
sociation resulting from the HL Tau UV field, although 
the uncertainties do not yet allow one to completely rule 
out standard ISM isotopologue ratios. Interpretation of 
these data are complicated by the fact that the line of 
sight toward HL Tau samples mostly the envelope rather 
than the disk. In addition, without the complementary 
data for C^^O, it is not possible to differentiate between 
mass-independent versus mass-dependent fractionation 
mechanisms; at temperatures that exist in the outer re- 
gions of disks, mass fractionation could be 20% or more 
(Young 2007), well within the range for HL Tau. 

Searching for oxygen isotope fractionation in proto- 
planetary disks requires an observational method that 
can measure the abundances of the four most common 
CO isotopologues to relative accuracies better than 10%. 
Isotopologue ratios measured using rotational emission 
lines depend on excitation and radiative transfer mod- 
els with signicant uncertainties due to non-thermalized 
level populations and differences in beam sizes arising 
from differences in the frequencies of the various CO 
isotopologues. Further, rotational emission is currently 
limited to the lowest rotational levels (J < 3), and is 
therefore tracing low-temperature gas only. An alterna- 
tive that eliminates most of these concerns is to use ab- 
sorption lines in which transitions from different isotopo- 
logues are known to probe nearly identical pencil-beam 
lines of sight. The fundamental ro-vibrational band at 
4.7/im and the overtone at 2.3 /xm of CO provide such 
an absorption-line tracer. We show that very high res- 
olution infrared spectroscopy can be used to derive ac- 
curate column densities of all 3 isotopologues of ^^CO. 
The main challenge is to locate sources with favorable ge- 
ometries providing column densities high enough to allow 
the detection of C^''0, but not so high as to extinguish 
the infrared source. One type of favorable absorption- 
line geometry is that of a disk viewed at high inclination 
angles to allow for absorption in the outer disk of in- 
frared continuum emission originating in the inner disk. 
In this paper, we take advantage of a new type of geom- 
etry, namely that of a binary star in which one star is 
being eclipsed by the disk of the other star. We compare 
one such binary T Tauri star with an embedded YSO 
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Fig. 1. — Schematic cartoon illustrating the two possible geometries for the VV CrA binary system. Case A illustrates the scenario 
whereby the disk of the secondary star is eclipsed by the outer disk of the primary; Case B shows the less probable, single inclined disk 
geometry. 



in which the absorption may in part be attributed to a 
remnant envelope. 



2. SOURCE SELECTION 

The first absorption-line source is the northern compo- 
nent of the binary T Tauri star VV CrA (J2000: RA = 
19 03 06.8, dec — —37 12 54.6); the southern component 
does not show CO in absorption. VV CrA N is part of 
the class of "infrared companion" (IRC) sources, which 
are highly extincted companions to optically visible T 
Tauri stars with projected separations of a few hundred 
AU or less (Koresko et al. 1997). Infrared companions 
typically show silicate absorption features, revealing that 
the absorbing material is localized. One interpretation is 
that the infrared companions have very compact, dusty 
envelopes. Here, we suggest a possible alternative expla- 
nation: that at least some of the infrared companions are 
regular T Tauri stars that are being eclipsed by the outer 
disk of the primary star of the binary (e.g., Hogerheijde 
et al. 1997). The only strict requirement for this scenario 
is that the disks of the two companions be non-coplanar, 
a geometry that gains support from the fact that non- 
coplanarity is a known property of some binary T Tauri 
stars (e.g., Koresko 1998; Stapelfeldt et al. 2003). The 
geometrical situation we envision is depicted schemati- 
cally in Figure 1. While the available data do not al- 
low a distinction between the two-disk scenario (Case A) 



and one in which the gas is located in the outer parts 
of a single inclined disk (Case B), the former scenario 
is geometrically more likely. In either case, we are ob- 
serving gas in the outer disk, and thus our analysis of 
the oxygen isotope ratios via CO is not predicated on 
the differentiation between the two geometries, although 
we note that Case B will minimize systematic differences 
between lines of sight at 2.3 and 4.7 /im. 

The projected separation of the VV CrA components 
is 2.1", corresponding to 273 AU at the distance of 130 
pc to the Corona Australis star-forming region (Casey 
et al. 1998; Neuhiiuser & Forbrich 2008). The VV CrA 
system is highly variable, as indicated by a K-band flux 
of -2.3 mag in 1987 (CheUi et al. 1995), 0.3 - 0.0 mag 
in 1993/95 (Ageorges et al. 1997; Chez et al. 1997) and 
1.3 mag at the time of our observation in 2007. Due 
to the fact that the southern component has remained 
relatively constant at K ^ 6 mag, this corresponds to the 
IRC fading by ^ 3 magnitudes in K in the past 20 years. 
Given the geometry we propose in Figure 1 (Case A), this 
variability could be a result either of intrinsic variability 
of the continuum source, or to column density variations 
in the orbiting disk of the companion that move across 
our line of sight to the continuum source. It is unlikely to 
be due to temperature or composition changes on such 
short time scales. 

The second source of this study is Reipurth 50 [RE 
50] (IRAS 05380-0728) (J2000: RA = 05 40 27.7, dec 
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TABLE 1 
Journal of Observations. 



Source/Obs. Date 


Spect. Range 


Int. time^ 




(UT) 


(^m) 


(minutes) 


S/Nb 


W Cr A/2007-08-31 


4.701 - 4.815 


32 


260 


W CrA/2007-08-31 


2.338 - 2.396 


23 


210 


Reipurth 50/2007-10-11 


4.645 - 4.901 


10-20 


290 


Rcipurth 50/2007-10-17 


2.338 - 2.397 


20 


40 



''Integration is variable with wavelength, so values are approxi- 
mate. 
''Median values. 



= -07 27 28), a luminous (250 Lq), embedded YSO in 
the Orion star-forming cloud at a distance of 470 pc. 
It illuminates a large, bright reflection nebula that ap- 
peared sometime between 1960 and 1970 (Reipurth & 
Bally 1986: Rcipurth & Aspin 1997). Strom & Strom 
(1993) found that Reipurth 50 is an FU Ori type YSO, 
and it has been shown to be surrounded by a compact 
0.1 Mq enyelope (Sandell & Weintraub 2001), placing it 
in stage I (Robitaille et al. 2006) of the YSO evolutionary 
sequence. It is likely that the absorbing gas is located in 
this envelope, rather than in a disk. Thus, Rcipurth 50 
and W CrA may provide a comparison between disk and 
envelope material, although larger samples of objects are 
required to make any firm conclusions. 

3. OBSERVATIONS 

The observations of W CrA and Reipurth 50 were 
taken as part of a large program to observe about 
100 YSOs and protoplanetary disks with the newly- 
implcmcntcd Cryogenic Infrared EchcUe Spectrograph 
(CRIRES) at the Very Large Telescope (VLT) in Chile. 
CRIRES is an AO-assisted spectrometer that operates at 
very high resolving powers (for infrared spectroscopy) of 
A/AA w 100000. The spectra of these sources allow for 
accurate determination of the CO column densities and 
isotopologue ratios. Table 1 summarizes the observing 
parameters for VV CrA and Rcipurth 50. 

The spectrum of VV CrA was observed on August 
31, 2007 in both the M (4.7 ^m) and K (2.3 ^im) bands. 
Reipurth 50 was observed on October 11 and 17, 2007, 
respectively, and the observations from both dates were 
combined. CRIRES M-band spectra of VV CrA obtained 
in April 2007 showed variations in the depths of the CO 
lines of ~ 10 percent relative to the August data, il- 
lustrating the importance of obtaining nearly concurrent 
K- and M-band spectra. The spectra of the v — {1 — 0) 
fundamental ro-vibrational bands for C^^O, C^®0, and 
i^CO and the v = {2-0) overtone band of C^^O for VV 
CrA were obtained using the 0.2" slit, resulting in a re- 
solving power of i? = A/AA « 95 000 (corresponding to 
3.16 km s~^), as measured on unresolved telluric lines. 
Several settings were obtained for each source to probe a 
range of rotational levels for all the isotopologues, span- 
ning J = to at least J = 9. For ^^CO, the v = (2-0) 
transitions were used because they have much smaller 
line strengths (by a factor of ~ 130) than the fundamen- 
tal, and are thus expected to be optically thin (or nearly 
so). 

The data were reduced using standard procedures, in- 
cluding flat-field correction, adjustments to account for 
detector non-linearity and linearization of the spectra in 
both the dispersion- and cross-dispersion directions. The 



spectra were wavelength calibrated using the telluric ab- 
sorption lines referenced to an atmospheric model spec- 
trum and transformed to the local standard of rest frame. 
Relative flux calibration was carried out by dividing the 
target sources by standard spectra of the stars HR 7236 
(B9V) and HR 1666 (A3III). For more information on the 
CRIRES data processing, see Pontoppidan et al. (2008). 

4. RESULTS AND DATA ANALYSIS 

The spectra of the fundamental and overtone CO 
bands are shown in Figures 2 and 3, respectively, il- 
lustrating the forest of narrow absorption lines due to 
CO gas. While we are only interested in the narrow 
CO absorption lines, it is noteworthy that VV CrA also 
shows broad, complex emission from CO, probably orig- 
inating in the innermost (i? < 1 AU) regions of the disk. 
Reipurth 50 docs not exhibit an emission component, but 
the known absorption band from solid CO centered on 
4.67y^m can be seen. Given that the respective emission 
and ice features in these sources are much broader than 
the absorption lines, their influence on derived gas CO 
column densities is minimal. In addition, the velocity 
shift of the center of the broad emission lines with re- 
spect to the narrow absorption lines are consistent with, 
and supportive of, our picture in which the absorption 
takes place in a different structure than the disk around 
the continuum source. 

Multiple lines from the four most abundant isotopo- 
logues of CO are detected in both sources, including the 
rare species C^^O (Figure 4). C^''0 has been observed 
in the ultraviolet toward the X Persei molecular cloud 
(Sheffer et al. 2002) and the P(l) C^^O ro-vibrational 
line has been recently reported in the T Tauri star SR21 
in infrared emission (Pontoppidan et al. 2008). While the 
latter observation would potentially allow a window into 
measuring isotope ratios in the inner <1 AU regions of 
disks and unambiguous determination of the location of 
the gas, radiative transfer modeling for emission lines is 
more complicated than for pencil-beam absorption lines. 
We emphasize that our observations of all four CO iso- 
topologues in infrared absorption probe a single sight 
line, therefore obviating the need for complex modeling 
required for emission lines. 

4.1. Derivation of column densities 

In order to determine the total column density for each 

isotopologue, the column density of CO in each rotational 
state must be measured; the total column density of CO 
can then be calculated by summing the measured col- 
umn densities in the observed ro-vibrational levels and 
assuming that the remainder are populated according to 
a Boltzmann distribution. If lines are spectrally unre- 
solved, accurate abundance derivations are complicated 
by the fact that it is difficult to determine optical depths, 
necessitating a curve-of-growth analysis. Given the high 
spectral resolution of CRIRES, the data reported in this 
study are not complicated by these factors, and we are 
able to distinguish intrinsic line broadening from that of 
the instrument. Figure 4 is a representative section of 
the fundamental CO band in Reipurth 50, illustrating 
the quality of the data. The CO lines were confirmed 
to be spectrally resolved by comparing them to adjacent 
unresolved telluric lines; the total CO line widths are 
slightly, but significantly, broader than the instrumental 
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Fig. 2. — Infrared absorption spectra of the CO fundamental ro-vibrational bands toward VV CrA and Reipurth 50. The fully observed 
spectral range is represented for VV CrA; for Reipurth 50, the portion with absorption lines used for the analysis is shown. Absorption 
lines are due to various CO isotopologues, representatives of which are indicated. Note that, in addition to the narrow absorption lines, 
VV CrA shows complex, low-level, broad-lined emission from hot CO gas, presumably located in the inner disk. The broad ^^C^^O ice 
feature in Reipurth 50 is also marked. 



resolution of 3.2kms~ exhibited by the telluric lines. 
Cohimn densities were obtained from the optical depths 
at line center (tq) derived from the best-fit profiles to 
each line. Tables 2 and 3 show the optical depths at line 
center and Doppler shifts (Vlsr) for each observed line 
used in the derivation of isotopologue column densities 
for VV CrA and Reipurth 50, respectively. 

The fundamental = (1 — 0) C^^O lines are optically 
thick. Therefore, the C^ O column densities were mea- 
sured using u = (2 — 0) transitions in order to use only 
lines that are as optically thin as possible. All line fits 
are consistent in having a single, constant width, within 
error, as expected for optically thin lines. That width is 
substantially larger than the thermal linewidth (see be- 
low), and so is likely due to some combination of local 
turbulent motions and radial velocity gradients trans- 



verse to the absorbing column. 

We fit the frequency-dependent flux, l^,^ for any given 
transition using 



ly = Ic exp(-s (j)^Nyj), (1) 

where Ic is the flux of the continuum and is the intrin- 
sic line profile (normalized to unity). The intrinsic line 
profile is assumed to be Gaussian with a full width at 
half maximum (FWHM) — 7intrinsic- N^j is the column 
density of the lower level of the transition, having rota- 
tional and vibrational quantum numbers J and v (v = 
for all observations presented here) , s is the integrated 
cross section, and To = sN^j(j)o- 
The integrated cross section (cm^s^^) (uncorrected for 
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Fig. 3. — Infrared absorption spectra of the CO overtone ro-vibrational band toward VV CrA [P(16)-P(17) not shown] and Reipurth 50. 
Detected ^'^C^'^O lines are marked. 




Wavelength [[im] 

Fig. 4. — Representative magnification of CO ro-vibrational lines 
for rare isotopologues in Reipurth 50. 



stimulated emission) is 



— ^ f 
rrieC 



(2) 



(Spitzer 1978), where fv'j'^vj is the oscihator strength 
for absorption from level u, J to level w', J', and is related 
to the Einstein A coefficient by 



A. 



v' J' -^vj 



9.J' 



fv'J' 



-vj 



(3) 



Einstein As were taken from the HITRAN database 
(Rothman 2005), which uses the values of Goorvitch & 
Chackerian (1994a, b). Ayres et al. (2006) compared dif- 
ferent sets of molecular parameters for the CO fundamen- 
tal band and found our adopted values to be consistent 
with other works. 

Because the intrinsic line widths are close to the in- 
strumental resolution, the Ii, given by (1) is convolved 
with the instrumental line profile modeled by a Gaus- 
sian with FWHM = 3.2 km s""'^ prior to fitting the line 
spectrum. We adopted the median intrinsic width of the 
C^^O lines as the best estimate of the line width and 
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Fig. 5. — Representative fits of CO isotopologues for VV CrA and Reipurth 50. Fits to C^''0 lines are to the overtone band. All other 
fits are to the fundamental spectra. The fits are shown here normalized to the linear best-fit baseline. Error bars are Iff/velocity channel. 
Statistical errors do not include telluric residuals (see text). 
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TABLE 2 

Optic;al depths and Doppler shifts for observed CO 

ISOTOPOLOGUE RO- VIBRATIONAL LINES USED IN DERIVING COLUMN 
DENSITIES IN VV CrA. UNCERTAINTIES ARE Icr DERIVED FROM 
THE LINE FITS. 
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()ptir;,l d.'prll Dn],pl.T ^.liifl 



Isot>,polnK,u. 


Li 




ID 




( 




) 








km 




12(^16o 


(2 





R(3) 





475 


± 





005 


4 


16 


± 


03 


12cl6o 


(2 





R{2) 





535 


± 





006 


4 


77 


± 


02 


12^160 


(2 










581 


± 





006 


5 


08 


± 


02 


12cl6o 


(2 


0) 


B(0) 





371 


± 





005 


5 


61 


± 


03 


12^160 


(2 


0) 


P(l) 





309 


± 





005 


5 


41 


± 


04 


12cl6o 


(2 


0) 


P(2) 





397 


± 





004 


5 


01 


± 


03 


12cl6o 


(2 


0) 


P(3) 





337 


± 





005 


4 


21 


± 


03 


12^160 


{2 


0) 


P(i) 





319 


± 





005 


3 


57 


± 


04 


12cl6o 


(2 





P(5) 





371 


± 





005 


3 


55 


± 


04 




(2 


0) 


P(6) 





381 


± 





004 


3 


67 


± 


03 


12cl6o 


(2 


0) 


P(7) 





360 


± 





005 


3 


62 


± 


03 




(2 


0) 


P(8) 





310 


± 





005 


3 


53 


± 


04 


12cl6o 


(2 





P(9) 





287 


± 





005 


3 


43 


± 


04 


12cl6o 


(2 





P(10) 





255 


± 





009 


3 


46 


± 


08 


12^160 


(2 





P(ll) 





231 


± 





005 


3 


78 


± 


05 


12cl6o 


(2 





P(12) 





206 


± 





004 


4 


08 


± 


05 


12cl6o 


(2 


0) 


P(13) 





154 


± 





005 


3 


99 


± 


08 


12cl6o 


(2 


0) 


f (14) 





136 


± 





005 


4 


78 


± 


09 


12cl6o 


(2 


0) 


P(15) 





088 


± 





007 


3 


94 


± 


16 


12cl6o 


{2 





P(16) 





093 


± 





007 


3 


72 


± 


17 


12cl6o 


(2 





P(17) 





064 


± 





005 


4 


28 


± 


21 


13cl6o 


(1 





H(7) 


1 


038 


± 





010 


3 


17 


± 


02 


13cl6o 


(1 


0) 


B(6) 


1 


381 


± 





006 


3 


25 


± 


01 


13cl6o 


(1 


0) 


H(5) 





880 


± 





007 


3 


08 


± 


01 


13cl6o 


(1 


0) 


K(4) 


1 


159 


± 





007 


3 


38 


± 


01 


13cl6o 


(1 


0) 


H(3) 


1 


104 


± 





008 


3 


30 


± 


01 


13cl6o 


{1 





P(2) 


1 


450 


± 





008 


4 


00 


± 


01 


13cl6o 


(1 





P(0) 





921 


± 





005 


5 


09 


± 


01 


13cl6o 


(1 





P(l) 





611 


± 





006 


4 


23 


± 


02 


13cl6o 


(1 


0) 


P(2) 


1 


025 


± 





010 


4 


40 


± 


02 


13cl6o 


(1 


0) 


P(4) 


1 


091 


± 





007 


3 


23 


± 


01 


12^180 


(1 





R(9) 





124 


± 





003 


2 


75 


± 


06 


12cl8o 


(1 





R(7) 





180 


± 





004 


3 


32 


± 


07 


12^180 


{1 





B(6) 





112 


± 





004 


2 


79 


± 


OS 


12cl8o 


(1 


0) 


K(5) 





179 


± 





004 


3 


38 


± 


05 




(1 


0) 


K(4) 





186 


± 





003 


3 


39 


± 


04 


12^180 


(1 


0) 


R(3) 





164 


± 





006 


3 


20 


± 


09 


12cl8o 


(1 


0) 


H(2) 





136 


± 





003 


3 


81 


± 


06 


12^180 


(1 










180 


± 





005 


5 


02 


± 


06 


12cl8o 


(1 





fl(0) 





141 


± 





004 




79 


± 


07 


12^180 


{1 


0) 


P(l) 





102 


± 





004 


5 


21 


± 


09 


12cl8o 


(1 


0) 


P(2) 





103 


± 





004 


3 


53 


± 


09 


12^180 


(1 


0) 


P(3) 





133 


± 





006 


3 


20 


± 


11 


12cl7o 


(1 


0) 


H(0) 





028 


± 





004 


5 


72 


± 


31 


12cl7o 


(1 





P(l) 





027 


± 





003 


5 


72 


± 


26 


12cl7o 


(1 





P(2) 





039 


± 





004 


4 


72 


± 


25 


12cl7o 


(1 





P(3) 





024 


± 





004 


3 


86 


± 


40 


12cl7o 


(1 


0) 


P(5) 





045 


± 





003 


2 


82 


± 


18 


12cl7o 


(1 


0) 


P(6) 





034 


± 





005 


3 


83 


± 


36 


12cl7o 


(1 


0) 


P(8) 





022 


± 





003 


3 


10 


± 


38 


12cl7o 


(1 





P(9) 





034 


± 





004 


4 


27 


± 


32 



fixed this parameter in the fits to individual lines. The 
measured mean FWHMs of the absorption profiles are 
4.5kms~-^ for Reipurth 50 and 5.1kms~^ for VV CrA, 
which, when deconvolved with the instrumental profile, 
yield 7intrinsic = 3.3 and 4.0 kms^-'^, respectively. The 
scatter of measured line widths suggests that the un- 
certainty on the median width is 10-20%. The value of 
7intrinsic IS rclatcd to thc Dopplcr velocity dispersion, b, 
by 7intrinsic = 2Vln2 b. The baseline intensities, Ic, were 
fit simultaneously with the line profiles using first-order 
polynomial fits to the local continuum. Sample fits are 
illustrated in Figure 5. 

From these fits we obtained rotational diagrams for 
the respective sources, shown in Figure 6. These plots 
show thc sub-lcvcl column densities, Nj/{2J + I), de- 
rived from the fits to the absorption lines, versus the 
energy of the state J. A linear trend is expected for 
a single-temperature gas. It is clear from the different 
slopes between the low- J transitions {J < 3) and higher- 
J transitions that there is a range of different gas tem- 



TABLE 3 

Optical depths and Doppler shifts for observed CO 

ISOTOPOLOGUE RO-VIBR ATIONAL LINES USED IN DERIVING COLUMN 

densities in Reipurtii 50. Un<:;ertainties are la derived 

FROM the line FITS. 



Reipurth 50 



IsPtnpnl.i^lIf 


Li 


ir 


ID 




Optir;,! -l.'pth 
I 'c.j 


Di 


ppl-'i 
.|, . k 


u 


llifl 


12^16o 


(2 





R 


3) 


1.013 ± 0.05 


5 


31 ± 





09 


12^160 


(2 





R 


2) 


1.035 ± 0.04 


4 


91 ± 





OS 


12(,16o 


(2 





R 


1) 


0.847 ± 0.04 


4 


94 ± 





10 


12cl6o 


(2 





R 


0) 


0.611 ± 0.03 


4 


35 ± 





11 


12^160 


(2 





p 


1) 


.619 ± .04 


4 


85 ± 





12 


12cl6o 


(2 





p 


2) 


0.775 ± 0.05 


4 


69 ± 





13 


12cl6o 


(2 


0) 


p 


3) 


0.799 ± 0.04 


5 


OS ± 





09 


12(-,16q 


(2 





p 


4) 


0.796 ± 0.04 


5 


36 ± 





10 




(2 





p 


5) 


0.752 ± 0.05 


4 


89 ± 





13 


12^160 


(2 





p 


6) 


0.658 ± 0.03 


4 


90 ± 





11 


12cl6o 


(2 


0) 


p 


7) 


0.609 ± 0.03 


5 


67 ± 





12 


12^160 


(2 





p 


8) 


0.486 ± 0.04 


5 


79 ± 





15 


12^160 


(2 





p 


9) 


0.291 ± 0.03 


5 


81 ± 





23 


12^160 


(2 





p 


11) 


0.220 ± 0.03 


5 


26 ± 





30 


12(-,16q 


(2 





p 


12) 


0.189 ± 0.03 


5 


58 ± 





33 




(1 





R 


15) 


0.314 ± 0.004 


5 


59 ± 





03 


13cl6o 


(1 





R 


13) 


0.558 ± 0.006 


5 


59 ± 





02 


13cl6o 


(1 





R 


12) 


0.746 ± 0.010 


5 


66 ± 





03 


13cl6o 


(1 


0) 


R 


11) 


0.881 ± 0.009 


5 


84 ± 





02 


13cl6o 


(1 





R 


10) 


1.044 ± 0.008 


5 


63 ± 





01 


13cl6o 


(1 





R 


7) 


1.451 ± 0.012 


5 


43 ± 





01 


13cl6o 


(1 





R 


6) 


1.510 ± 0.006 


5 


44 ± 





01 


13cl6o 


(1 





R 


5) 


1.877 ± 0.007 


5 


13 ± 





00 


13cl6o 


(1 





R 


4) 


1.955 ± 0.007 


5 


26 ± 





00 


13cl6o 


(1 





R 


3) 


2.099 ± 0.008 


5 


26 ± 





00 


13cl6o 


(1 





R 


2) 


2.449 ± 0.011 


4 


94 ± 





01 


13cl6o 


(1 





R 


0) 


1.840 ± 0.006 


4 


74 ± 





00 


13cl6o 


(1 





P 


1) 


1.612 ± 0.005 


5 


17 ± 





00 


13cl6o 


(1 





P 


4) 


1.978 ± 0.012 


5 


08 ± 





01 


13cl6o 


(1 


0) 


P 


9) 


1.244 ± 0.008 


3 


61 ± 





01 


13cl6o 


(1 


0) 


P 


10) 


1.087 ± 0.008 


5 


56 ± 





01 


13cl6o 


(1 





P 


11) 


0.897 ± 0.009 


6 


09 ± 





02 


13cl6o 


(1 





P 


13) 


0.611 ± 0.007 


4 


35 ± 





02 


12cl8o 


(1 





R 


15) 


0.025 ± 0.003 


5 


56 ± 





30 


12<-,18o 


(1 





R 


14) 


0.051 ± 0.004 


5 


67 ± 





19 


12cl8o 


(1 





R 


13) 


0.084 ± 0.006 


6 


02 ± 





IS 


12(,18o 


(1 





R 


12) 


0.101 ± 0.005 


6 


05 ± 





12 


12(^180 


(1 





R 


11) 


0.142 ± 0.004 


5 


66 ± 





06 


12(-.18q 


(1 





R 


10) 


0.196 ± 0.003 


5 


44 ± 





04 




(1 


n 


R 


9) 


0.239 ± 0.004 


5 


40 ± 





04 


12(,18o 


(1 





R 


7) 


0.372 ± 0.003 


5 


55 ± 





02 


12(-,18o 


(1 





R 


6) 


0.450 ± 0.003 


5 


52 ± 





01 


12cl»0 


(1 





R 


5) 


0.485 ± 0.003 


5 


78 ± 





01 


12(-.18q 


(1 





R 


4) 


0.457 ± 0.003 


5 


21 ± 





02 


12clSo 


(1 





R 


2) 


0.651 ± 0.003 


5 


08 ± 





01 


12^180 


(1 





R 


1) 


0.568 ± 0.004 


5 


13 ± 





01 


12cl8o 


(1 





R 


0) 


0.348 ± 0.002 


5 


15 ± 





01 


12^180 


(1 





P 


2) 


0.456 ± 0.004 


5 


15 ± 





02 


12cl8o 


(1 





P 


3) 


0.516 ± 0.005 


5 


21 ± 





02 


12pl8o 


(1 





P 


8) 


0.312 ± 0.004 


3 


90 ± 





03 


12cl8o 


(1 





P 


9) 


0.238 ± 0.006 


5 


57 ± 





06 


12^180 


(1 





P 


10) 


0.193 ± 0.004 




98 ± 





05 


12cl7o 


(1 





R 


9) 


0.049 ± 0.003 


5 


36 ± 





15 


12^170 


(1 





R 


8) 


0.068 ± 0.003 


5 


29 ± 





10 


12cl7o 


(1 





R 


7) 


0.082 ± 0.004 


5 


68 ± 





11 


12cl7o 


(1 





R 


6) 


0.097 ± 0.006 


5 


78 ± 





13 


i^ci^o 


(1 





R 


5) 


0.114 ± 0.006 


5 


87 ± 





11 


12cl7o 


(1 





R 


4) 


0.128 ± 0.004 


5 


27 ± 





07 


12cl7o 


(1 





R 


2) 


0.132 ± 0.003 


4 


80 ± 





05 


12cl7o 


(1 





R 


0) 


0.082 ± 0.003 


5 


10 ± 





OS 


12^170 


(1 





P 


2) 


0.115 ± 0.003 


5 


OS ± 





06 


12cl7o 


(1 





P 


3) 


0.105 ± 0.002 


5 


25 ± 





05 


12c"o 


(1 





P 


4) 


0.097 ± 0.003 


5 


36 ± 





06 


12cl7o 


(1 





P 


5) 


0.106 ± 0.003 


5 


48 ± 





06 


12cl7o 


(1 





P 


8) 


0.054 ± 0.005 


6 


07 ± 





19 


12cl7o 


(1 





P 


9) 


0.041 ± 0.003 


5 


47 ± 





17 



peratures along the line of sight toward both sources. 

This is to be expected for essentially all geometries in 
which the gas is heated by the central star. While as- 
suming a continuous distribution of temperatures might 
be most realistic, we found that a distribution consist- 
ing of only two temperatures is sufficient to capture the 
distribution of points in the excitation diagrams for both 
sources. In the absence of any explicit information on the 
actual temperature distributions, we proceeded with a si- 
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TABLE 4 

Measured CO isotopologue rotational temperatures.'' 





VV CrA - high T 


VV CrA - low T 


Reipurth 50 - high T 


Reipurth 50 - low T 


12^160 


261 ±5 


16 ± 1 


152 ±9 


19 ±3 


13^160 


> 100 


16 ±2 


212 ±4 


18 ± 1 


12C1H) 


176 ± 11 


6 ± 1 


150 ± 3 


21 ± 2 




17:-! ± 30 


l.j 


110 (> 


Ki ± 3 



^Values were obtained from a simultaneous two-temperature 
model. All temperatures in K. 



TABLE 5 

Measured CO isotopologl'e c;olumn densities and ratios. Uncertainties reflect the 68% confidence level 



VV CrA - high T VV CrA - low Reipurth 50 - high T Reipurth 50 - low T'' Local ISM'' 
N(i2ci6o) 7.2±0.2 X 10i*cm-2 1.4 ± 0.1 x lO^* cm-^ 9.0 ± 0.5 x lO^* cm-^ 1.6 ± 0.2 x lO^* cm-^ 
N(i2ci6o)/N{i3ci6o) lOOdilO'^ 110 ±8 110 ±7 65 ± 7 69 ± 6 

N("c"'0)/N(i2ci**0) 690 ± 30 1700 ± 250 490 ± 30 340 ± 40 557 ± 30 

N(i2ci6o)/N(i2ci70) 2800 ±300 4200 ± 800 2200 ± 150 2300 ± 340 2005 ± 155 

N(^^Ci^O)/N(i^C^^O) 4.1 ±0.4 2.5 ±0.6 4.4 ± 0.2 6.7 ± 0.8 3.6 ± 0.2 



''Low-temperature values are shown for completeness; see text for 
discussion on their uncertainties. 
''(Wilson 1999) 

"assuming TC^'^CO) = T(^^CO) 

multaneous two-temperature fit to the data for each iso- 
topologue (Figure 6); from these models, individual low 
and high temperatures were derived. For both objects, 
the two temperatures consisted of a cold component of 
T~ 10-20 K and a warm component of T ~ 150-250K. 
In the case of VV CrA, it is worth noting that the tem- 
perature of the warm component is roughly consistent 
with that of the uppermost surface layers at a few 100 
AU in detailed disk heating models (Kamp & DuUemond 
2004; Jonkheid et al. 2004). The column density of the 
warm component is significantly larger than expected 
from these models perpendicular to the midplane, but 
may be explained if small dust grains have been depleted 
from the upper disk layers with a concommitant drop in 
opacity. A non- vertical, slanting line of sight through the 
disk can also enhance the total column of warm gas. Fur- 
ther modeling of these specific sources is needed to deter- 
mine if the observed temperatures and column densities 
are quantitatively consistent with flared disk models. 

In both VV CrA and Reipurth 50, the derived low 
temperatures are consistent with the outer disk near the 
midplane or, in the case of Reipurth 50, possibly the sur- 
rounding envelope. The high-temperature contributions 
to the measured absorption strengths of lines dominated 
by the cold component are considerable (Figure 6). 

The uncertainties in individual line depths are domi- 
nated by systematics, and these are taken into account in 
the final errors for the isotopologue ratios. The greatest 
source of systematic uncertainty is likely to be residu- 
als from telluric lines that were not completely removed. 
This source of error is difficult to quantify, and varies 
strongly for difi^erent lines. To account for this, we added 
5% uncertainties in quadrature to each line strength for 
most of the isotopologues in both temperature regimes; 
exceptions to the 5% uncertainty were accounted for in 
VV CrA for the low temperature ^^C^^O, low temper- 
ature ^^C^^O and high temperature ^^C^^O lines, for 
which no systematic uncertainty, 10 % and 20 % sys- 
tematic uncertainty were included, respectively. These 
variations resulted from variable optimizations in fits to 



obtain reduced chi-square values near unity; chi-squared 

ellipses arc shown in the tcmpcrature-column density 
plane (Figures 7 and 8). A few lines are clearly out- 
liers relative to a linear fit in the rotational plot by many 
standard deviations; these lines, shown as faded points in 
the fits in Figure 6, were excluded from the final fits. In- 
clusion of the estimates for systematic errors means the 
uncertainties cited here are likely to be conservative. 

4.2. Isotopic Ratios and Uncertainties 

Derived excitation temperatures and isotopologue ra- 
tios are listed in Tables 4 and 5, respectively. There is 
a significant difference in excitation temperatures among 
the different isotopologues as evidenced by the variable 
slopes in Figure 6 for each temperature component, ap- 
parently violating the assumption of local thermody- 
namic equilibrium. For example, the excitation tem- 
perature for the high-temperature ^^CO component in 
Reipurth 50 is greater than that of the other isotopo- 
logues (Figure 8). A similar effect is seen in the warm 
VV CrA ^^CO lines (Figure 7). However, the appar- 
ently higher temperature of ^^CO could be due to line 
trapping in optically thick ro- vibrational lines, in which 
energy escapes from the rarer isotopologues more readily 
than from ^^CO due to high optical depth in the latter. 
In this case, the derived isotopologue ratios should not 
be affected. 

Because the ^^CO lines are observed in a different 

wavelength region (2.3 /im) than the other isotopologues 
(4.7/xm), it is possible that the source location and size 
are not exactly the same due to the greater scattering 
of the K band relative to the M band. This could cause 
significant differences in the column densities that are 
probed. If the absorbing gas is very close to an in- 
frared emitting source that changes in size or position 
with wavelength, it is still possible that wc are probing 
different volumes and perhaps path lengths. However, 
in a spherical envelope, material with temperatures of 
150-200 K is located at distances of -100-200 AU from 
a 250 Lq source, compared to a size of only a few AU for 
the infrared source. Hence, this is potentially the largest 
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Fig. 6. — Rotational excitation diagrams for CO isotopologues 
in VV CrA and Reipurth 50. Simultaneous fits to the two- 
temperature model and derived high and low temperatures for each 
isotopologue are indicated by the red curves and dashed Hnes, re- 
spectively. Error bars are la propagated from the Gaussian fits. 
Ej is the energy of the Jth rotational state above the ground rota- 
tional state and k is the Boltzmann constant. Faded symbols have 
been excluded from the fits. 

source of error for both an embedded source with a large 
reflection nebula like Reipurth 50, as well as a single in- 
clined disk geometry as illustrated by Case B (Figure 1) 
for VV CrA. The more probable two-disk scenario for VV 
CrA (Case A) is not likely to suffer from optical path dif- 
ferences. With CRIRES, the pointing is the same for all 
spectral settings, with both K- and M-band observations 
tracking sources in the K band. We therefore believe 
that there are no significant differences in the absorbing 
volumes for the K and M bands for Reipurth 50 and VV 
CrA. 



The low temperatures derived from our model are low 
enough that they could well correspond to the tempera- 
ture of the surrounding molecular cloud for Reipurth 50, 
or the outer disk in Case A or B for VV CrA, perhaps 
weakly heated by the nearby protostar. However, due 
to the significant contribution of the high-temperature 
component on the low-temperature lines, we believe that 
our low-temperature ratios contain a much larger degree 
of uncertainty than we are able to precisely quantify. We 
therefore do not believe these numbers are meaningful. 
In discussing isotope ratios, we therefore refer to the ra- 
tios obtained from the high-temperature lines only. 

5. DISCUSSION 

Both VV CrA and Reipurth 50 should have 
[Ci60]/[Ci80] and [Ci60]/[Ci^0] values similar to the 
local ISM in the absence of extensive isotope partition- 
ing by photochemistry (self-shielding) or mass-dependent 
kinetics. The oxygen isotope ratios for the ISM as a 
function of distance from the Galactic center (Galacto- 
centric radius, Rqc) can be estimated from the 4 kpc 
and 8 kpc values reported by Wilson (1999), yielding for 
[i60]/[iS0]: 

[i60]/[i«0]j,,^^ = (57.5 ± 10) Rgg + 97, (4) 

where the la uncertainty in the slope reflects the uncer- 
tainty in the means (standard errors) reported in Table 4 
of (Wilson 1999). For [i60]/[i'^0] we can combine equa- 
tion 4 with values for [i^0]/["0] in the ISM. These es- 
timates vary from 3.5 (Wilson 1999) to 4.1 (Wouterloot 
et al. 2008). 

Based on equation 4 and [i*0]/[^'^0] = 4.1 for 
the ISM, one should expect CO surrounding VV CrA 
(i?c;c = 7.8kpc) to have [Ci^OJ/pi^O] of 550 ± 90 and 
[Ci60]/[Ci^0] of 2300 ± 250. The measured values of 
690±30(lcr) and 2800±300 for VV CrA are greater. Con- 
versely, the expected [Ci^Ol/pi^O] and [Ci^OJ/pi^Q] 
values for Reipurth 50 {Rgg = 8.4kpc), 580 ± 90 (Icr) 
and 2400 ± 270, respectively, are indistinguishable from 
the measured values of 490 ± 30 (Icr) and 2200 ± 150 
(Table 5). 

Results for the oxygen isotopologues in the high- 
temperature regime are summarized in the three-isotope 
plot shown in Figure 9, where differences in [i^O]/[i^O] 
and [^^0]/[^^0] from local ISM values are shown in per 
mil using the linearized form of the delta notation com- 
monly used in isotope cosmochemistry. Here, 6^^0' — 
103ln(i8i?,/i8i?LocaiiSM) and iSiJ. = [Ci8o]/[ci60]. 
Error ellipses in the three-isotope plot were calculated 
using a Monte Carlo calculation that propagates un- 
certainties in column densities to the final delta val- 
ues. The three-isotope plot shows that the differences 
in [if50]/[i80] and [^'^0]/["0] between best estimates 
of the local ISM and the disk surrounding VV CrA are 
not consistent with mass-dependent isotope fractionation 
(represented by a line with a slope of 0.52 in Figure 9). 
This conclusion is firm to the extent that it can be as- 
sumed that the precursor molecular cloud to VV CrA 
was similar to the local ISM. The conclusion is robust 
despite the uncertainty in the [i^0]/[^^0] ratio of the 
ISM because the position of the VV CrA data point on 
the higher end of the range in reported [i'*0]/[i^O] val- 
ues (near 4.1) precludes mass-dependent fractionation as 
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Fig. 7. — Goodness-of-fit contours for the high temperature lines (J > 4) of the four measured isotopologues in VV CrA. The contours 
show the 68, 95, and 99% confidence levels. Red stars indicate total column densities obtained from the simultaneous two-temperature 
model. 



a primary control on the CO isotopologue ratios in all 
cases. Rather, the trajectories between VV CrA and 
plausible local ISM isotopic compositions have slopes of 
1 or greater. The slope-1 relationship between local ISM 
and CO in the VV CrA disk is reminiscent of the mass- 
independent fractionation that characterizes oxygen in 
the solar system and is a tell-tale sign of self-shielding 
during CO photolysis. 

The ability to distinguish slope-1 (mass-independent) 
from slope-0.52 (mass-dependent) trends in Figure 9 is 
critical for interpretation of the data and is a key ad- 
vantage of having data for all three oxygen isotopologues 
of CO. The deficit of the rarer isotopologues, C^^O and 
C^^O, relative to C^^O in the VV CrA disk, is on the 
order of 20 to 40% respectively, which is significant at 
the > 2(7 level. 

The [i2C0]/["C0] ratios for both VV CrA and 
Reipurth 50 are similar, with high-temperature values of 
100 ± 10 and 110 ± 7, respectively (Table 5). Ratios for 
both objects are high compared with the local ISM value 
of 69±6 (Wilson 1999). In a case where CO self-shielding 
was almost certainly responsible for extreme C^^O ex- 
cesses in the study of the X Persei molecular cloud by 
Sheffer et al. (2002), there are no associated excesses in 
^^C, suggesting that carbon and oxygen isotope effects 
during CO photolysis can be decoupled. Because both 
of our objects have the same [^^C]/[^^C] but one has an 
excess in C^^O relative to the ISM, our data also suggest 
a decoupling of C and O isotope ratios in CO in these 
objects. Self-shielding by CO during photolysis should 
in principle result in an excess of ^^CO relative to ^'^CO. 



However, atomic carbon liberated by CO photolysis is 
photoionized largely to C+, and 0+ reacts rapidly to ex- 
change carbon with CO by the reaction ^^0+ -|- ^^CO ^ 
12Q+ _^ 13QQ (warin et al. 1996; van Dishoeck & Black 
1988; Langer et al. 1984), hence driving the ^^C back 
into CO and diminishing the signature of selective pho- 
todissociation. In contrast, oxygen liberated by CO pho- 
tolysis remains neutral in the photodissociation regions 
(PDRs) and is relatively unreactive; oxygen will not ex- 
change with CO or other molecules prior to sequestration 
as H2O produced either in the gas phase by a reaction 
channel involving H;^ (Herbst 2000) or on grain surfaces 
(Hasegawa et al. 1992). The decoupling between carbon 
and oxygen effects seen here can be explained in the con- 
text of photochemical self-shielding but offers as yet no 
explanation for the anomalously high [^^CO]/[^^CO] for 
both VV CrA and Reipurth 50. This anomaly will be 
the focus of a later publication. 

The cause of the disparities in [^^0]/[^®0] and 
ji6Qj^ji7Qj between the disk surrounding VV CrA and 
the local ISM cannot yet be deduced unequivocally, as 
more data points for Figure 9 are needed. However, the 
mass-independent character of the displacement of this 
source is most consistent with a photochemical enhance- 
ment of C^^O relative to both C^^O and C"0 due to 
self-shielding by the more abundant C^^O isotopologue. 

The absence of a similar effect in Reipurth 50 might 
be attributable to the different geometry of the absorb- 
ing component of this source. The protostellar enve- 
lope in Reipurth 50 is presumably less affected by CO 
self-shielding than the disk in VV CrA. If the difference 
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Fig. 8. — Goodness-of-fit contours for the high-temperature lines (J > 4) of the four measured isotopologues in Reipurth 50. The contours 
show the 68, 95, and 99% confidence levels. Red stars indicate total column densities obtained from the simultaneous two-temperature 
model. 
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Fig. 9. — Comparison of oxygen isotope ratios between the in- 
terstellar medium and CO surrounding the two YSOs reported 
here (stars). Ratios for the YSOs are derived from the high- 
temperature CO abundances. The solid and dashed lines show 
mass-independent fractionation lines with different assumptions for 
the ^^O/^^O ratio, as indicated. The grey line shows the mass- 
dependent fractionation line. Ellipses represent 95% confidence 
limits derived from the goodness-of-fit contours shown in Figures 
7 and 8. 

in mass-independent isotope fractionation between disks 
and protostellar envelopes indeed holds for more objects, 
an important question is what is the cause. The time- 
scale for the isotope-selective photodissociation and sub- 



sequent chemistry itself should be short, less than 1000 
years, in the PDR of both types of objects. However, the 
lines of sight likely probe more than just the PDR layer 
so the isotope signature of these other layers is equally 
relevant. In the disks, the photodissociation is confined 
to thin surface regions but the ^^0-depleted atomic oxy- 
gen can be mixed vertically down to the midplane region 
on timescales of '^lO^ years in the outer disk, where it 
can then be sequestered into H2O (Young 2007). To- 
gether with the effects of a non-vertical geometry (§4.1) 
and grain growth by coagulation, this can result in an 
isotopic fractionation effect over a much larger column 
than would be expected from a thin PDR layer; grain 
growth can lead to higher UV flux through the disk and 
greater photodissociation. In contrast, Reipurth 50, be- 
ing only in stage I of its evolution, may not have had 
sufficient time to acquire a measurable C^^O excess due 
to mixing (~ 10"" yr), so that, in the absence of mixing, 
the isotope fractionation of the PDR layer contributes 
little to the total column. 

6. CONCLUSIONS 

Our new high-resolution IR absorption data for two 
YSO environments, a protoplanetary disk and a proto- 
stellar envelope, demonstrate that measurement of both 
[Ci60]/[Ci80] and [Ci60]/[Ci^0] in CO is possible with 
sufficient precision to distinguish photochemical effects 
from mass-dependent isotope fractionation. The new 
high-resolution data exhibit a detectable deficit in C^^O 
and C^^O in the disk surrounding VV CrA relative to the 
local interstellar medium. The CO surrounding Reipurth 
50 exhibits no discernible differences in [C^^O]/[C^*0] 
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and [Ci^O]/[Ci^O] relative to the local ISM. One likely 
explanation for the different oxygen isotopologue ratios 
for these objects is that a photochemical deficit in C^^O 
and C^^O relative to C^^O proceeds only in a disk ge- 
ometry and may require > 10^ years in the environments 
surrounding YSOs. Future analyses of more objects in 
different stages of evolution are now needed to enable dif- 
ferentiation between potential self-shielding in the disk 
versus the parent cloud. 
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